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Abstract 

Neutrino mass sum rules are an important class of predictions in flavour models relat¬ 
ing the Majorana phases to the neutrino masses. This leads, for instance, to enormous 
restrictions on the effective mass as probed in experiments on neutrinoless double beta 
decay. While up to now these sum rules have in practically all cases been taken to hold 
exactly, we will go here beyond that. After a discussion of the types of corrections that 
could possibly appear and elucidating on the theory behind neutrino mass sum rules, we 
estimate and explicitly compute the impact of radiative corrections, as these appear in 
general and thus hold for whole groups of models. We discuss all neutrino mass sum 
rules currently present in the literature, which together have realisations in more than 
50 explicit neutrino flavour models. We find that, while the effect of the renormalisation 
group running can be visible, the qualitative features do not change. This finding strongly 
backs up the solidity of the predictions derived in the literature, and it thus marks a very 
important step in deriving testable and reliable predictions from neutrino flavour models. 
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1 Introduction 


The last two decades have lifted neutrino physics to new heights, as experiments have consid¬ 
erably driven the field. Nowadays it can be seen as an established fact that the neutrino mass 
and flavour bases are different 1 . In a nutshell, a certain flavour (say, an electron neutrino 
z/ e ) does not have a well-defined mass but is instead a superposition of three active neutrino 
mass eigenstates ^ 1 , 2 , 3 - Mathematically, this change of the basis is described by a mixing 
matrix called the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix, 


t^PMNS 


R 23 U 13 R 12 P 0 

( C12C13 
— S12C23 — Cl2S23Sl3e 1<5 
S12S23 — Ci2C23Sl3e 1<5 


S12C13 

C12C23 — Sl2S23'Sl3e 1 ' 5 
—C12S23 — •Sl2C23Sl3e 1<5 


Si3e~ l5 \ 

S23C13 Po, 
C23C13 J 


( 1 . 1 ) 


where <5 is the Dirac CP-phase and Pb = diag(e -1<?;>1//2 , e - 1 < ^ 2 / 2 ,1) is a diagonal matrix containing 
the two Majorana phases <^ 1 , 2 - While experimentally we have a relatively clear picture about 
mixing in the lepton sector, in the sense that we have by now measured all three mixing 
angles $12,13,23 to some precision, we have no idea why their values are so large compared to 
those of the quark sector [2 . The most popular idea to explain these obvious patterns is to 
relate them to the properties of discrete symmetry groups [3|, although alternative ideas such 
as, e.g., a random mass pattern 4j or a transmission from other sectors 5 do exist. 

One basic problem of neutrino flavour models based on discrete symmetries is their in- 
distinguishability at low energies: would the models predict mixing angles far from their 
experimental values, we would consider the models to be excluded. If, however, they predict 
values within the current 3a ranges, we have a hard time distinguishing them unless the 
experimental precision is considerably improved, which is not to be expected very soon. If 
on the other hand the models predict correlations between certain low-energy observables, 
these could be used as additional handles. Well-known examples for testable correlations 
are neutrino mixing sum rules [6j, such as S 23 — ^ = — ^-cosd, which, e.g., allow to pre¬ 
dict the Dirac CP phase from some of the mixing angles. The type of correlation we would 
like to study here, instead, are so-called neutrino mass sum rules , which relate the three 
(complex) neutrino mass eigenvalues to each other. Mass sum rules have been studied for 
quite some time (7-12|, but it was only in recent years that systematic analyses have been 


presented [13 15]. These analyses clearly show that, among all observables, the effective neu¬ 
trino mass |m ee | as measured in neutrinoless double beta decay (0 is(3/3) 116] can be modified 
most significant!}^] - in a way that several groups of models could be excluded in the near 
future, despite the uncertainties involved, in particular if in addition to a new limit or even 
a measurement of 0 vf3f3 information on the neutrino mass ordering (i.e., whether normal, 
mi < m 2 < m 3 , or inverted, m 3 < m\ < m 2 ) was available [15]. Turning the logic round, 
future experiments could even “gauge” their definitions of stages with increasing exposure 
using the predictions from sum rules |17|. 


Mhe reason for other observables like the effective electron-neutrino mass square or the sum E of all 
light neutrino masses not to be affected very strongly is that they contain no Majorana phases and thus are 
only constrained by “half” of the information contained in a neutrino mass sum rule, which is intrinsically 
a complex equation. Furthermore, both these observables are comparatively insensitive to changes in the 
smallest neutrino mass, if that is small by itself; they may however be affected in cases where a sum rule 
forbids a certain mass ordering. 
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What all previous studies have in common is that they treated neutrino mass sum rules 
as if they were exact to each order, i.e., as if they would be perfectly known0 This leads 
to relatively strong predictions such as some sum rules excluding a particular mass ordering 
even in the case where neutrinos are very close in mass, such that the two orderings should 
be hardly distinguishable 115]. However, in general neutrino mass sum rules are very unlikely 
to hold exactly, since several types of corrections could appear. In particular two types of 
corrections are evident: 


• First of all, mass matrices in concrete models are typically only computed at leading 
order. Thus, when taking into account next-to-leading order (NLO) corrections, a 
mass sum rule may be destroyed (although counterexamples are known [l8j). Further 
corrections may arise, too, e.g. from normalising non-canonical kinetic terms correctly. 
However, all these types of corrections are difficult to analyse in a unified manner, since 
they are very model-dependent. Furthermore, it is not a priori clear how to treat the 
case of an “approximate” sum rule in an accurate way. 


Second, as any quantity in a quantum field theory, the values of neutrino masses vary 
depending on the energy scale, a fact known as renormalisation group evolution (RGE) 
or, in a less formal manner, simply dubbed running. The running may also affect the 
validity of a sum rule, and in particular it may modify the allowed regions and/or 
open up or close down the consistency of the sum rule for a particular mass ordering. 
The good point with these types of corrections is that, while also they are in principle 
model dependent, one can at least study their effect on classes of models which can be 
effectively described by the Weinberg operator 19 , i.e., where all right-handed neutrinos 
are so heavy as to be integrated out at a relatively high energy scale, such that their 
exact mass spectrum does not affect the low-energy mass matrices significantly. Note 
that we assume neutrinos here to be Majorana particles. 


In this work we will focus on the effect of the second type of corrections, the radiative cor¬ 
rections, since for them practical consequences can be derived for many cases without having 
to resort to a model-by-model investigation. While this approach of course implies that we 
intrinsically disregard the first type of corrections, we would like to stress once more that 
in some cases the NLO corrections do not modify the sum rules |18| . While up to now no 
general criteria for this behaviour to appear are known, it is at least clear that it can poten¬ 
tially happen in which case the RGE corrections are the only general corrections that exist. 
Thus we will in what follows investigate the effect that RGE corrections have on the ranges 
predicted for \m ee \ from neutrino mass sum rules. 

This paper is organised as follows: in Sec. [2] we shortly review how to derive predictions 
from neutrino mass sum rules, before explaining in more detail both the general effect of 
radiative corrections and our numerical approach in Sec. |3j The results for all known sum 
rules, along with an accompanying discussion, can be found in Sec. [4j We then conclude in 
Sec. [5} Some subtleties related to computing roots of complex numbers, which are decisive 
in order to derive the correct predictions from those sum rules involving square roots, are 
discussed in Appendix |A| 

2 Ref. 13 did study the effect of possible deviations of a few per cent, however, in that case the deviation 
was assumed to be proportional to one of the light neutrino masses, which is not only unmotivated but also 
introduces an undesired measure-dependence. 
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2 Reviewing neutrino mass sum rules 


In this section we want to briefly review how neutrino mass sum rules can be parametrised 
and how they can be interpreted as a prediction for the Majorana phases as a function of 
the (physical) neutrino masses. A very detailed description can be found in Ref. |l5| , while 
a more pedagogical introduction to the broader topic of neutrino flavour models, featuring 
several example models that partially also lead to sum rules, can be gained by studying the 
known reviews (3j. 

The essential feature of a flavour model incorporating a (light) neutrino mass sum rule is 
that the eigenvalues of the light neutrino mass matrix depend on two complex parameters only. 
Typically this can arise from any neutrino mass generation mechanism in which the structure 
of one mass matrix is generated by two flavor^ couplings while all other matrices only have a 
single scale which can be factored outj^] For illustrational purposes we will discuss now exactly 
such a two flavon case, where we suppose that in a certain model, the light neutrino mass 
matrix m u is generated by a type I seesaw mechanism 21 , such that m v = 


l R m D 

with m,D {Mr) being the Dirac (right-handed neutrino) mass matrix. While in the most 
general case, the Dirac mass matrix mp is practically arbitrary (apart from being at most of 
electroweak size) and the Majorana mass matrix Mr is only symmetric, in a concrete flavour 
model their structure may be further constrained. For example, the different generations of 
right-handed neutrinos and charged leptons could be chosen to transform under particular 
representations of a family symmetry, such that the mass matrices are given by 


/ 1 

0 

0 \ 1 

f 2 a s + do 

cx s 

-cts \ 

m D = 0 

0 

1 yv and Mr = 

-a. 

2a s 

—a s + do 1 

V o 

1 

0 / 

V ~a s 

—a s + ao 

2 a s 1 


A 


( 2 . 1 ) 


where y is a Yukawa coupling, v is the electroweak VEV, and A is a the mass scale of the right- 
handed neutrinos (see Ref. [l5| for details). Note that Mr does have a particular structure 
depending on the two dimensionless couplings d Sj o- These couplings arise as ratios of two 
different flavon VEVs and the breaking scale of the family symmetry, and the representation 
of the flavons determines in which entries of Mr the two parameters show up. The Dirac mass 
matrix mj}, in turn, does not involve any flavon coupling and its size is entirely determined 
by the product of the Yukawa coupling and the electroweak VEV. While it does have some 
non-trivial structure owing to the family symmetry, the mass scale yv can be factored out of 
this matrix. This is the decisive point: looking at the resulting light neutrino mass matrix, 


m v = 


y 2 v 2 


a 0 (a 0 + 3d s )A 


( cto + ot s 
a s 
a s 


Q. s 

d s (l - 3b) 
a s + boiQ 


a s \ 
a s + ha o I , 
a s (l - 3b)) 


( 2 . 2 ) 


with b = ao/(cto — 3d s ), depends only on two paramaters (namely a s and do) in what its 
structure is concerned, while any mass scales can be factored outj^] Computing the complex 

3 A flavon is a scalar field that is a total singlet under the standard model gauge symmetry, but it trans¬ 
forms non-trivially under the family symmetry and thus breaks it spontaneously when obtaining a vacuum 
expectation value (VEV). 

4 We want to remark that there are also other cases imagineable and indeed known. For instance, the 
model 20] has only one flavon relevant for the light neutrino masses which receives a VEV which depends on 
two complex mass parameters. 

5 Such a structure can only be achieved if exactly one type of matrix in the seesaw formula is generated by 
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eigenvalues (mi, m 2 , m 3 ) = ^ 3ct 1 <77 > 3q , j of this neutrino mass matrix, one can 
see immediately that a sum rule of the form 


1 

rhi 


1 

777,3 


2 

777-2 


(2.3) 


is implied (see Sec. 4.7 for a discussion of exactly this sum rule). Note that, in fact, the 
only true requirement on the light neutrino mass matrix to fulfill a neutrino mass sum rule 
is the dependence on exactly two parameters, up to an overall scale. While of course this 
structure is somewhat related to the mixing pattern predicted by a certain model, there is 
no direct constraint that the mixing would have on the sum rule. In particular, part of the 


mixing could be induced by the charged lepton sector which is not involved in Eq. (2.2). 


Note futher that, at least in the case at hand, there is in fact also a mass sum rule for the 
eigenvalues (Mi, M^, M 3 ) = (ceo + 3a s ,o:o, —ao + 3a s )A of the right-handed Majorana mass 
matrix Mr: Mi — M 3 = 2 M 2 , again because it depends on exactly two parameters up to an 
absolute scale. While such a relation might have interesting implications at high energies, it 
is experimentally hardly accessible at low-energies, as long as the right-handed neutrinos are 
very heavy. As a final note, we should mention that typically the parameters playing the role 
of a s and ao are complex numbers, which is exactly what yields to a neutrino mass sum rule 
imposing a relation between the two light neutrino Majorana CP phases, while however their 


exact values are not constrained since the prefactor in Eq. (|2.2|) will be complex in general. 
The most general mass sum rule can be written as0 

(2.4) 


Ai 777 ,fe 1X1 + A 277 rfe 1X2 + A 3 ?ri 3 e 1X3 = 0 


where rhi labels the complex mass eigenvalues (i.e., including the Majorana phases), d 7 ^ 0, 
and Ai > 0. The phase Xi £ [0,27t) originates from the sum rule itself, i.e., it contains both 
the phase of Ai and a possible minus sign. Following similar steps as in 


15 


we can rewrite 

Eq. (2.4) in a more convenient form. First we express the complex masses in terms of the 
Majorana phases 07 G [0, 2n) and the physical mass eigenvalues mi > 0 &s rhi = 777 je 1Qi . 
Dividing Eq. (2.4) by A 3 and using the abbreviations c- L = Ail A 3 > 0 and ai = Xi + & a i, we 
end up with: 


ClTTlfe 1 " 1 + C27772e 1Q:2 + Ttlfe'" 3 = 0 . 
Now we can multiply Eq. (2.5) by e _1Q3 = e~ 1 (x 3 +da 3 ) and define 

= Xi~ X3 ■ 


(2.5) 


( 2 . 6 ) 


With the use of the Majorana phases —cbi = d* — d 3 , which appear in Eq. (1.1), we end up 
with our final sum rule: 


S(777l, 777 2 ,777 3 , 01,02) D, C 2 , d, A\'i 3 , Ay' 23 ) = 

ci ^777ie _1< ^ e lAxi3 + C 2 (m 2 e -102 ) i 


,iAx23 


+ 777 3 = 0 . 


(2.7) 


precisely two flavon couplings, while all other mass matrices need at most one flavour (so that the VEV can 
be factored out). 

6 Note that we are using conventions different from the ones used in 15 . We do this in order to match the 
conventions used in REAP/MPT 22 , which is our tool of choice to compute the RGE corrections. 
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^‘ 01 +A ^13 


Figure 1: Geometrical interpretation of the mass sum rule. 


Similarly as in 13-15 , we interpret the mass sum rules geometrically, since this equation 


describes the sum of three vectors which form a triangle in the complex plane, see Fig. [lj Via 
the law of cosines we can express the angles a = —dcf) 2 + A Y 23 and f3 = —dcj) 1 + Ayi 3 in terms 
of the masses: 


cos a = 


cos f3 = 


c\m\ d — — m\ d 

drmd 


2 c 2 m?,m 3 


— c\rri\ d — m\ d 
dZ 2 , 


2c-imf m 3 


( 2 . 8 ) 

(2.9) 


These equations decide about the validity of a mass sum rule since the right-hand side has 
to be in the interval [— 1 , 1 ] to obtain real values for a and /3 or, respectively, 4 > 1 and cj) 2 . 
Since the cosine is an even function we obtain two solutions for by taking its inverse. 
This is connected to the fact that the orientation of the triangle in the complex plane is 
irrelevant. Nevertheless one encounters a subtlety when the sum rule involves square roots of 
the masses. Since the square root of a complex number is not uniquely defined, the orientation 
of the triangle in the complex plane in fact turns out to be important (see Appendix [A|. 

The basic reasons for the rise of sum rules is that, in most neutrino mass models, the 
structure of the physical light neutrino mass matrix arises from products of several mass 
matrices. The most generic example is again a type I seesaw mechanism |21 , where the light 
neutrino mass matrix m v = is generated from a multiplication of the Dirac 

mass matrix mo and the heavy neutrino mass matrix Mr- If, in this product, the structure 
of mo is generated by two flavon couplings, whereas the scale of Mr can be factored out, 
this would lead to a sum rules featuring a square root, such as y/rh\ ± \Jfh 3 = 2 \fni2- On 
the contrary, if the structure of Mr is generated by two flavon couplings, whereas mo only 
features a single scale, the resulting sum rule would feature an inverse power of one, like 
2 / 777.2 = 1/m-i + l/m 3 or l/rh\ + l/rfi 2 = I/ 7713 . This scheme extends to basically all kinds 
of neutrino mass models, see Ref. 115] and in the last Ref. |3| for more detailed explanations. 
In Tab. □ we have collected all the sum rules we found in the literature and which we will 
discuss in the following with their parameters c\ , C 2 , d, A\i 3 , and A \'23 that characterise them 


according to Eq. (2.7). 
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Sum rule 

References 


Cl 

C2 

d 

A\'13 

AX23 

1 



23 

-29 


1 

1 

1 

7r 

7r 

2 



30 



1 

2 

1 

7r 

7r 

3 | 

1C 

1,24 

-28 

5,31 

33 

1 

2 

1 

7r 

0 

4 



34 



1/2 

1/2 

1 

7T 

7T 

5 



35 



2 

x/3+1 

y/3-l 

a/3+1 

1 

0 

7r 

6 

tJI 

,18 

20 

23 

36 

1 

1 

-1 

7r 

7r 

7 

|8- 

11 

32, 

33, 

37 

1 

2 

-1 

7r 

0 

8 



38 



1 

2 

-1 

0 

7T 

9 



39 



1 

2 

-1 

7T 

7t/2,37t/2 

10 


1 

2,4 

0 


1 

2 

1/2 

7T, 0, 7t/2 

0,7r, 7 t/2 

11 



14 



1/3 

1 

1/2 

7r 

0 

12 



41 



1/2 

1/2 

-1/2 

7r 

7r 


Table 1: Summary table of the sum rule we will analyse in the following. The parameters 
Ci, C 2 , d , Ayi 3 , and A \23 that characterise them are defined in Eq. (2.7). In sum rule 9 and 10 
two possible signs appear which lead to two possible values of A.Xi 3 - Note that sum rule 10 
with Ayi 3 = A\23 = 7t/2 is the sum rule in 


12 


which was wrongly interpreted before. 


We do not want to explicitly discuss any models which lead to the sum rules here. Note 
that, however, there is no model which predicts sum rule 11, but it was shown in 14 that 
such a sum rule can be predicted in models with a type I seesaw mechanism. 


3 The implications of renormalisation group running and how 
to compute them 

In this section, we will estimate how renormalisation group running may affect a neutrino mass 
sum rule. In particular, we will answer the question whether it is possible for a radiatively 
corrected sum rule to allow for mass orderings which are forbidden at tree-level (i.e., if the 
exact sum rule holds exactly). We then discuss why for some sum rules we can expect sizeable 
RGE effects even for a small mass scale which one would not expect from the RGE running 
of the parameters itself. But before we come to these two points, we start with some general 
remarks. 

3.1 The general effect of radiative corrections 

The running of neutrino masses and mixing parameters is already known for quite some time, 
see, e.g., (42|. Naively one might expect that the running of the mixing parameters is small 
and that visible effects only happen if we have a large Yukawa coupling. In the SM there is 
no large Yukawa coupling in the lepton sector but in the MSSM, for large tan /3, the Yukawa 
couplings can even be of 0( 1). This expectation will be confirmed later where we find only 
small corrections in the SM, while for the MSSM with tan/3 > 20 the corrections start to 
become interesting. 

On top of this effect there can be an additional enhancement of the RGE effects induced 
by a parametric enhancement. Some corrections, e.g., for the mixing angles and phases, are 
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proportional to m 2 /Am 2 , where m labels the lightest neutrino mass and Am 2 is one of the 
neutrino mass squared differences. In the quasi-degenerate mass regime this easily yields an 
enhancement of 0 ( 100 ). 

As discussed in |42|, the masses themselves run mostly due to the Higgs wave function 
renormalisation which includes the top Yukawa coupling but which is flavourblind and there¬ 
fore will not matter for us. But the tan /3 enhancement and the parametric enhancement for a 
large neutrino mass scale will directly induce large RGE effects for the Majorana phases, such 
that their low-energy values can be very different from the constrained high-energy values Q 

This can be easily seen in our numerical results later on but before we get there we want to 
discuss two other questions which can be understood better by estimates instead of extensive 
numerical parameter scans. 


3.2 Trying to reconstitute forbidden mass orderings 

Some sum rules are only viable for one mass ordering at tree-level as it was already noted 
before 13-15]. For instance, in what we labelled sum rule 2, where (ci,c 2 ,d, A\i 3 , Ay 23 ) = 
(1, 2,1, 7 r, 7 r), only normal mass ordering is allowed. A natural question is whether RGE 
corrections can reconstitute the inverted ordering. To answer this question, we can have a 


closer look at Eq. (2.8). For the RGE-corrected value of cos a, we expand the masses and 
find: 


cos cr ree + <5(cosa) = 


(mi + dmi) 2 — 4(?b 2 + dm 2) 2 — (m 3 + dm 3 y 
4 (m 2 + 5m. 2 )(m 3 + dm 3 ) 


2 /t ™2 _ (m 2 + 4 m 2 — m|) 


m 2 — 4 m 2 — m 3 mi 

4m 2 m,3 2?n 2 m3 1 

(mi 2 — 4m 2 2 + m 3 2 ) 


4m 2 ?B3 


dm2 


4m 2 m| 


-dm 3 


(3.1) 


where the rrq are the low energy neutrino masses and the dmi their respective RGE corrections. 
We want to obtain inverted ordering, where m 3 < m\ < m 2 . This implies that 


m 2 — 4m| — m 3 < — (3m| + m 3 ) A 


1 


< 


1 


cosa tree = 


m 2 m-3 mg 
( m\ — 4m 2 — m|) 1 

4m 2 ?7i3 4 



(3.2) 

m 2 \ 

3 —| + 1 ) < — 1 • 

. rn 3 J 

(3.3) 


Thus, inverted ordering is ruled out on tree-level or low energies. Now, if d(cos a) is sufficiently 
positive at high energies, we might nevertheless realise this regime. The correction dmi to 

7 There is some subtle issue involving the mass ordering. Just by the running of the masses itself the 
ordering will hardly flip. Nevertheless, by running 612 , for instance, could turn negative at high energies which 
one might compensate by exchanging the first and second mass state and hence A?n|i becomes negative at 
the high scale which would also lead to kinks and jumps in the running of the angles and phases. To avoid 
confusion with this effect we have chosen conventions where the mass hierarchy is always preserved. 
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the z-th neutrino mass eigenvalue can be estimated as 42 


1 


5mi = 16vr2 (“ RGE + 2C yr S 12 s h) m l lo g , 

1 («RGE + 2 Cylc\ 2 slz) m 2 log + 0(6 13 ) , 


5m 2 = 
5m 3 = 


167T 2 

1 

167T 2 


(«rge + 2Cylc \ 3 ) m 3 log + 0 ( 6 1 3 ) , 


(3.4) 

(3.5) 

(3.6) 


where we have assumed the parameters to be constant at leading order, so that we integrate 
the j3 function between the Z-scale Mz and [i > Mz- Note that «rge ~ 3 is a function of 
gauge and Yukawa couplings, while C = — 3/2 in the Standard Model (SM) and C = 1 in the 
minimal supersymmetric Standard Model (MSSM). 

We can plug this back into Eq. (3.1) and use 0 23 « 7 t/4 and sin 6*12 ~ l/\/3 to find: 


<5(cosa) ~ — 


Cy^ 3m 2 — 4to 2 + m§ // 


192vr 2 


rn 2 m :i 


3 log —— 
6 M z 


(3.7) 


The first thing to note is that the dependence on orge drops out (which is true for all sum 
rules). Because of 3m 2 —Am^+m^ < 0, the corrections are negative for the MSSM and hence 
they make cos a even smaller. For the SM, on the other hand, they would have the right sign 
- but then we would need ig ^ 2 log to be of 0(1), which implies ^ to be way beyond the 
Planck scale. Hence, we conclude that the inverted ordering cannot be reconstituted by RGE 
corrections for the second sum rule. 

Similar estimates can be done for the other sum rules with missing mass orderings (i.e., 
sum rules 2, 3, 4, 5, 10, 12). In fact, for sum rule 12 (ci = 1/2, C 2 = 1/2, d = —1/2, 
Ayi 3 = 7 r, A \'23 = 7r), the MSSM corrections would have the right sign so that we could hope 
to reconstitute the missing ordering in that case, but according to our estimate we would 
need for a neutrino mass scale below 1 eV a tan j3 value of more than 500, which practically 
excludes this possibility as well. 


3.3 Impact of the RGE corrections for a small mass scale 


By looking at the formulas for the RGE effects on the masses and hence on cos cc, one might 
think that they have barely an impact for a small mass scale since the RGEs are proportional 
to the mass scale itself. But we will show that indeed also a small mass scale can lead to 
significant corrections for cos a, which happens due to a parametric enhancement of the form 
\JAm 2 /m 2 that is large for small masses. 

As an example we have a look at sum rule 1 with an inverted ordering (m = m 3 ). The 
expression for cos a for sum rule 1 reads 


cos cr ree + 5 (cos a) 


2 2 2 
nil — m 2 — m 3 

2 m 2 m 3 


+ 


Ai log jL 

96-7T 2 Mz 


—3 m\ + to 2 — mjj 


m 2 m 3 


(3.8) 


where we have already used the previous estimates and d 23 « n /4 and sin6\ 2 ~ l/\/3. 
Expressing the masses in terms of m 3 via the mass squared differences and neglecting small 
terms of order m^/Am 2 , we end up with a negative tree-level value: 


cos a 


tree 


Am \ 

2m 3 y / | Am| 2 


(3.9) 
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For the correction term we find: 


<5(cos a) 


C Vr 1q Jf_ ( 2 \ Am h\ \ 
967 T 2 & Mz y m 3 -^/| Am 2 2 1 J 


(3.10) 


From the tree-level term we get a lower bound on m 3 , given by m 3 = 7.6 • 10 -4 eV. The 
correction has a negative sign in the MSSM - just as the leading order term - which means 
that the corrections increase the lower bound on the masses for this sum rule. For the mass 
where cos a tree = — 1 , the correction further decreases the value of cos a and hence it will be 
forbidden. In addition, we see that the corrections are enhanced for a small mass scale due 
to the small values in the denominator. Therefore the total cos a gets very sensitive to small 
changes in m 3 (due to RG corrections) and hence the allowed range for the Majorana phases 
at the high scale becomes larger. 

This correction differs from the leading order term by a factor of A = 4 log ' 

To get a better feeling for the size of the effect of the corrections we plug in tan f3 = 50, 
m = 8 • 10 -4 eV, and fj, = Ms = 10 13 GeV to get: 


A « 0.83 . (3.11) 

This implies an increase of the lower bound of m 3 by 83%. For sum rule 4 we find in a similar 
way a 83% correction to m 3 , as well. Other sum rules in both orderings do not cover such 
small mass scales, and hence we do not find such an enhancement of the RGE corrections for 
them. 

Later on, in our numerical scans, we will obtain a lower bound for m 3 for sum rule 1 of 
about m ~ 9.1 • 10 ~ 4 eV, which is a much smaller effect than our estimate suggests. But 
indeed for the parameter point we mentioned our one-step approximation for the /3 functions 
is not very good because also the angles, especially 612 , run significantly. Nevertheless, with 
our estimates we can easily understand the apparent broadening of the allowed region for sum 
rules 1 and 4 which we will see later on. 


3.4 The numerical approach 


Since the constraints on the mixing parameters are satisfied at different energy scales, their 
experimental ranges, cf. Tab. [2j restrict the mixing angles and the mass squared differences 
at a low energy scale M^j^] 

In contrast, the mass sum rules in fact constrain the Majorana phases as functions of the 
lightest neutrino mass which we will label in the following simply as m at a high energy scale, 
where we assume the sum rule to be predicted by the respective flavour model. As a generic 

10 13 GeV, and we then employ a 


choice we set this scale equal to the seesaw scale 
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running procedure between Ms and Mz- Choosing a different scale instead would not change 
our results dramatically, as long as it is not different from Ms by several orders of magnitude, 
so that the running would extend over a considerably larger or smaller energy range. 

In our scans we will present results for the SM extended minimally by the Weinberg oper¬ 
ator to accommodate neutrino masses, as well as for the MSSM plus the Weinberg operator. 


s To be precise, the actual experiments performed detect neutrinos of even lower energies, between MeV and 
GeV depending on the source. However, it is a well-known fact that the change of these parameters between 
the scale Mz and low energies is negligible, provided that the particle content is that of the SM. 
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Parameter 



best-fit (Tier) 

3 <7 range 

0\2 in 

O 




33-481°;?® 

31.29 -> 35.91 

CO 

5’ 

O 




8.501°!? © 8.511°!? 

7.85 —>9.10 ©7.87 —^ 9.11 

623 in 

O 




42.31?;° © 49 . 5 I 2 ! 

38.2 —»• 53.3 © 38.6 —»• 53.3 

5 in 0 





251 +67 
/0i -59 

0 -> 360 

A m 2 2l 

in 

10 -5 

eV 2 


7 5f) +0 ' 19 
' ' OU -0.17 

7.02 -> 8.09 

Amli 

in 

10" 3 

eV 2 

(NO) 

2 4 57+ 0 - 047 

Z.UOf-0 047 

2.317 —»• 2.607 

Am .32 

in 

10 -3 

eV 2 

(IO) 

—2 449+ 0 - 048 

-2.590 -> -2.307 


Table 2: The best-fit values and the 3 a ranges for the parameters taken from [l]. The two 
minima for both $13 and 6*23 correspond to normal and inverted mass ordering, respectively. 


The most relevant supersymmetry (SUSY) parameter for the running is tan/3, which we have 
chosen to be 30 or 50 in our scans, while the exact mass spectrum of the SUSY particles plays 
hardly any role. We have fixed the SUSY scale, where we switch from SM to MSSM RGEs, 
to 1 TeV but - again - the dependence on this scale is only logarithmic and hence very weak. 
Furthermore we have neglected the SUSY threshold corrections for the masses and mixing 
parameters j43|. Both for small tan/3 and in the SM the running is small, and hence the 
results in these two cases would be very similar. In fact the SM results will look very similar 
to the results without RGE effects at all, cf. Ref. 15 , as to be expected from the small size 
of the relevant corrections. 

In order to perform our numerical computations, we have made use of the REAP/MPT 
package 1221. To do this we run the parameters up to the high scale Mg and calculate there 


the modulus of the sum rule [i.e., of the left-hand side of Eq. (2.7)], which is minimised with 
respect to the low energy Majorana phases such that, if the sum rule holds, the minimum of 
the modulus should yield a numerical zero. The mixing angles and mass squared differences 
are varied within their experimental 3cr-rangesJ^] while 5 and m are free parameters at the low 
scale. We vary the Dirac CP phase 5 between 0 and 2ir, since it has not been measured yet 9 10 
and we have also scanned over values for the lightest mass between 1 • 10 -4 and 0.15 eV. The 
upper bound on m is chosen in accordance with the cosmological bound on the sum of the 
neutrino masses |44j: 

0.17 eV, (3.12) 


although we should note that what is displayed is the average limit taken between the two 
mass orderings. 

One might wonder if the running of the Majorana phases is sufficiently large to alter the 
results obtained in previous studies [l5|. We will show that the running of the Majorana 
phases is negligibly small in the SM, whereas we do see a substantial effect in the MSSM. 
Furthermore we will show that the running is not identical for (j)\ and <f> 2 - 


9 In principle we could have employed a similar procedure using their best-fit values at low energies. However, 
since the continuous comparison between low- and high-energy scales is numerically expensive, this would blow 
up the computational time without significant gain. 

10 Note that the global fits indicate some finite range at lcr level, however, this is not a significant tendency 
at the moment. 
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The running of the Majorana phases is given by 1221: 


; Cy T ( misf 2 sin + m 2 cf 2 sin c ^2 , mim 2 cf 2 s | 3 sin(0i - 0 2 ) 

01 = -jzjt ( m 3 cos(20 23 )- 1 ^-h 


02 = 


47r 2 

Cyr 

47r 2 


Am| 2 


Am|i 


\ m i s i 2 sin 01 + m 2 c? 2 sin 0 2 m 1 m 2 sf 2 s| 3 sin(0i - 0 2 ) 

m 3 cos(2w 23 )-—- x3 -—-h 


A m § 2 


Ama 


+ 0 ( 013 ) , 

(3.13) 

+ 0 ( 013 ) , 

(3.14) 


where we have used the abbreviations s 




= sin and c,; 
A m2 


uij = cos Vij, and we have neglected 
factors proportional to the small quantity ^'"p 1 • The two formulas are identical except for 

the factor c 2 2 for 0i instead of s\ 2 for 0 2 in the last term. This small difference is nevertheless 
crucial for the difference in the running of the phases. Since cf 2 is about two times larger 
than sf 2 for the best-fit value of 0 32 , and since this term is additionally enhanced compared 


to the first term in (3.14) due to the small mass difference in the denominator, the running 
of 0i is considerably stronger than the running of 0 2 . With increasing mass scale, the RGE 
effects drive 0i 2 to smaller values. This further enlarges the difference in the running of the 
Majorana phases, since sf 2 is decreasing whereas cf 2 is increasing with smaller values of 0i 2 . 

Due to the enhancement in the /3-functions of the MSSM governed by tan /3, we see this 
difference best in the running of the phases in the MSSM. As a typical example we compare 
the predicted low scale values of the Majorana phases as a function of m coming from sum 
rule 5 within the SM and the MSSM for tan /3 = 50 in Fig. [2] The black lines represent the 
predicted values for the phases without taking the RGE corrections into account, while the 
red points are the results from our numerical approach as described above. For 0i the red 
points strongly deviate from the black lines whereas the points for 0 2 gather around the black 
lines which supports our argument. 

Since the Majorana phases themselves are not directly measurable in the near future, 
we will present in the following section our results in terms of predictions for the allowed 
range of the effective neutrino mass \m ee \ as potentially measured in Qu/3/3. This observable 
is explicitly given by: 

\m ee \ = Imit/g 2 ! + m 2 t/ 2 2 + m 3 t/ 2 3 | = mic 2 2 c 2 3 e+ 1?: ’ 1 + m 2 s 2 2 c 2 3 e -19i2 + m 3 sl 3 e~ 2lS . (3.15) 

In all cases, as explained above, we will compute the predictions for the SM and for the MSSM 
(with tan/3 = 30 and 50), the latter of which can lead to considerably different predictions. 


4 Numerical results for concrete sum rules 

In this section we employ the procedure as described in the previous section to obtain allowed 
ranges for the smallest neutrino mass eigenvalue m and \m ee \ for all sum rules we found in 
the literature. Note that the numerical values obtained in this section may be limited by the 
finite statistics of our numerics. Furthermore, the oscillation parameters have been updated 
since the time when Ref. [ 15] has been written, which accounts for the small differences we 
obtain compared to that reference. In our plots regions with inverted mass ordering are drawn 
in yellow while regions with normal mass ordering are drawn in blue. 
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Figure 2: Predicted values of the Majorana phases in the SM (upper plots) and in the MSSM 
with tan/? = 50 (lower plots) as a function of m (which is m 3 in the case of sum rule 5). The 
black lines represent the predicted low scale values of the Majorana phases without taking 
RGE corrections into account, while the red points are the results of our numerical approach. 


4.1 Sum Rule 1: mi + m 2 = m 3 

The parameters for this sum rule are (d, c\, C 2 , Axi 3 , AX 23 ) = (1,1,1, vr, 7r), and the corre¬ 
sponding plots look like: 
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m [eV] 


This sum rule yields (ra m i n , |m ee | m in) = (0.026,0.026) eV ((0.00065,0.015) eV) for normal 
(inverted) mass ordering, if running with the SM particle content is applied, which is consistent 
with the values obtained in Ref. 115] (see discussion in Sec. 7.7 therein). For tan/3 = 30 (50), 
the values change to (0.028,0.026) eV ((0.028,0.026) eV) for NO and to (0.00065,0.014) eV 
((0.00075,0.015) eV) for 10, respectively. 


4.2 Sum Rule 2: ui! = m 3 — 2m 2 

The parameters for this sum rule are (d, ci, c 2 , Axi 3 , A% 23 ) = (1,1, 2 , n, ir), and the corre¬ 
sponding plots look like: 





Disfavoured by Qv/3{3 



This sum rule predicts normal ordering only, and with the SM particle content it yields 
(T7i min . |ra ee |min) = (0.016,0.015) eV, if the running is applied. These numbers are consistent 
with the values obtained in Ref. 115] (see discussion in Sec. 7.10 therein). For tan/3 = 30 
(50), the values basically remain at (0.016,0.015) eV ((0.016,0.015) eV), while still only NO 
is allowed. 


4.3 Sum Rule 3: mi = 2m 2 + m 3 

The parameters for this sum rule are (d, ci, c 2 , A%i 3 , Ay 23 ) = (1,1, 2,7r, 0), and the corre¬ 
sponding plots look like: 
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m [eV] 


This sum rule predicts normal ordering only, and with the SM particle content it yields 
(m min , |m ee |min) = (0.016,0.0042) eV, if the running is applied. These numbers are consistent 
with the values obtained in Ref. 15 (see discussion in Sec. 7.9 therein)]^] For tan/3 = 30 (50), 
the values change to (0.016,0.0036) eV ((0.016,0.0036) eV), while still only NO is allowed. 


4.4 Sum Rule 4: fri! + m 2 = 2m 3 

The parameters for this sum rule are (d, ci, c 2 , Axi 3 , A\' 2 3 ) = (1,1/2,1/2, n, 7r), and the 
corresponding plots look like: 



m [eV] m [eV] m [eV] 

This sum rule predicts inverted ordering only, and with the SM particle content it yields 
('/Timm-1?^ee|min) = (0.00028, 0.015) eV, if the running is applied. These numbers are consistent 
with the values obtained in Ref. [15] (see discussion in Sec. 7.12 therein). For tan/3 = 30 (50), 
the values change to (0.00030,0.014) eV ((0.00040,0.014) eV), while still only IO is allowed. 

4.5 Sum Rule 5: mi — ^~ 1 m 2 + x/ ^ +1 m 3 = 0 

The parameters for this sum rule are (d, ci, c 2 , A%i 3 , A% 23 ) = (1, and the 

corresponding plots look like: 

11 However, note that our computation just misses the cancellation region, in contrast to the one presented 
in Ref. 15 . Nevertheless there is no real discrepancy, since the question whether or not all parameters can 
conspire to yield \m ee \ practically zero depends strongly on the actual oscillation parameters used [45], and 
our values are updated compared to the ones uses in publications two years ago. 
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m [eV] 


This sum rule predicts inverted ordering only, and with the SM particle content it yields 
(rri m in , |m ee | m i n ) = (0.024,0.053) eV, if the running is applied. These numbers are consistent 
with the values obtained in Ref. [l5| (see discussion in Sec. 7.14 therein) f^j For tan (3 = 30 
(50), the values practically remain at (0.024,0.053) eV ((0.024,0.053) eV), while still only IO 
is allowed. 


4.6 Sum Rule 6: The sum rule 1/mi + l/m 2 = l/fri 3 

The parameters for this sum rule are (d, ci, C 2 , Ayi 3 , Ax' 23 ) = (—1,1,1, n, 7 r), and the corre¬ 
sponding plots look like: 



This sum rule yields (m m i n , |m ee |min) = (0.010,0.0016) eV ((0.028,0.048) eV) for normal 
(inverted) mass ordering, if running with the SM particle content is applied, which is consistent 
with the values obtained in Ref. 115] (see discussion in Sec. 7.1 therein). For tan/3 = 30 (50), 
the values change to (0.011,0.0017) eV ((0.011,0.0017) eV) for NO and to (0.028,0.052) eV 
((0.028,0.054) eV) for IO, respectively. 


4.7 Sum Rule 7: 1/mi — 2/m 2 — l/m 3 = 0 

The parameters for this sum rule are (d, ci, C 2 , A%i 3 , AX 23 ) = (—1,1, 2,7r, 0), and the corre¬ 
sponding plots look like: 

12 Note the typo in the value for m m i n in table 4 of that reference. 
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This sum rule yields (m m i n , |m ee | m in) = (0.0044,0.0046) eV ((0.017,0.019) eV) for normal 
(inverted) mass ordering, if running with the SM particle content is applied, which is consistent 
with the values obtained in Ref. 115] (see discussion in Sec. 7.8 therein). For tan/3 = 30 (50), 
the values change to (0.0044,0.0045) eV ((0.0044, 0.0047) eV) for NO and to (0.017,0.018) eV 
((0.017,0.017) eV) for 10, respectively. 


4.8 Sum Rule 8: 2/m 2 = I/ur + l/m 3 

The parameters for this sum rule are (d, ci, c 2 , A%i 3 , A% 23 ) = (— 1,1, 2 , 0, ir), and the corre¬ 
sponding plots look like: 



1 




Disfavoured by Qv/3{3 



This sum rule yields (m m i n , |m ee | m i n ) = (0.0044,0.0045) eV ((0.017,0.015) eV) for normal 
[inverted] mass ordering, if running with the SM particle content is applied, which is consistent 
with the values obtained in Ref. 115] (see discussion in Sec. 7.6 therein). For tan/3 = 30 (50), 
the values change to (0.0044,0.0044) eV ((0.0044, 0.0047) eV) for NO and to (0.017,0.019) eV 
((0.017,0.018) eV) for 10, respectively. 


4.9 Sum Rule 9: l/m 3 + 2i(—l) ?? /m 2 = 1/fhx 

The parameters for this sum rule are (d, ci, c 2 , A\i 3 , A\' 2 3 ) = (—1, l,2,7r,7r/2 or 37r/2), de¬ 
pending on whether rj = 0 or 1, and the corresponding plots look like: 
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For both rj = 0,1, this sum rule yields (m m i n , |m ee | m in) = (0.0044,0.0028) eV ((0.017, 0.016) eV) 
for normal (inverted) mass ordering, if running with the SM particle content is applied, which 
is consistent with the values obtained in Ref. {l5| (see discussion in Sec. 7.2 therein). For 
tan/3 = 30 (50), the values change to (0.0044,0.0028) eV ((0.0044,0.0030) eV) for NO and to 
(0.017,0.017) eV ((0.017,0.016) eV) for IO, respectively. 

This may at first look surprising, however, even though the RGEs for the Majorana phases 
(and hence the corresponding predictions) are different in both cases, this information gets 
lost when varying over the Dirac CP phase <5, as we have checked numerically. Turning the 
argument round, if 5 was known at least to some extend, we would potentially be able to 
distinguish the two versions of this sum rule. 


4.10 Sum Rule 10: \/m^ + i n ^/m 3 = 

The parameters for this sum rule are (d, ci, C 2 , Axi3, Ay 2 . 3) = (1/2,1, 2, it or 0 or n/2 ,0 or it or 7 t/2 ), 
depending on 77 = 0,1, 2, and the corresponding plots look like: 
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For each value of rj. this sum rule predicts normal ordering, and with the SM particle content it 
yields (m m i n , |m ee |min) = (0.00093,0.000014) eV, if the running is applied. These numbers are 
consistent with the values obtained in Ref. |15| (see discussion in Secs. 7.4 and 7.13 therein). 
For tan (3 = 30 (50), the values change to (0.00093,0.000025) eV ((0.00093,1.5 • 10 -6 ) eV), 
while still only NO is allowed. 

4.11 Sum Rule 11: 3^m 2 + 3\/m 3 = \J mi 

The parameters for this sum rule are (d, ci, c 2 , Axi 3 , A% 23 ) = (1/2,1/3,1,7r, 0), and the cor¬ 
responding plots look like: 



m [eV] 



m [eV] 


in [eV] 


This sum rule yields (m m i n , |m ee | m i n ) = (0.032,0.022) eV ((0.024,0.041) eV) for normal (in¬ 
verted) mass ordering, if running with the SM particle content is applied. Since this sum rule 
is hypothetical, in the sense that no explicit underlying model is known yet, no numerical 
predictions have been listed in Ref. 115] (see discussion in Sec. 7.5 therein). However, the 
left plot seems quasi identical to the one presented in that reference. For tan/3 = 30 (50), 
the values change to (0.032,0.021) eV ((0.032,0.021) eV) for NO and to (0.024,0.044) eV 
((0.024,0.042) eV) for IO, respectively. 



4.12 Sum Rule 12: 1 /= 2/ \/m^ — 1 / \Zm -2 

The parameters for this sum rule are (d, c±, c 2 , Axi 3 , A% 23 ) = (—1/2,1/2,1/2,7r, 7 t), and the 
corresponding plots look like: 
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This sum rule predicts normal ordering, and with the SM particle content it yields (m m i n , |m ee |min) 
(0.0027,0.0031) eV, if the running is applied. These numbers are consistent with the values 
obtained in Ref. 115] (see discussion in Sec. 7.3 therein). For tan/3 = 30 (50), the values 
change to (0.0027,0.0032) eV ((0.0027,0.0032) eV), while still only NO is allowed. 


4.13 Discussion 


As can be seen, RGEs can have a non-trivial effect on the regions allowed by certain sum rules. 
Although we “only” presented scatter plots (for a good reason though, see the discussion in 
Sec. 3.4), the tendency is clearly visible. In most cases (i.e., sum rules 1, 2, 4, 5, 6, 7, 8, 11, 12) 
the effect of the RGEs is to broaden the allowed regions, although in one case the broadening 
occured only within the parameter region that is already disfavoured by current neutrino 
mass bounds (sum rule 5) and in two cases it only appeared for inverted mass ordering, since 
only a small mass range for the normal ordering is allowed in these sum rules where we do 
not have an enhancement effect of the RGEs (sum rules 7 and 8). In a few cases (sum rules 
3, 9, 10), there is no visible effect in the plots, even though - numerically - the parameters 
do run. These three sum rules at least naively seem to have nothing in common, so that a 
simple “accidental” cancellation of the running effects is unlikely. Rather, there is a more 
fundamental reason: in these sum rules, large RGE effects are suppressed by the values of the 
angles and the phases at the high scale. 


Furthermore, as anticipated in Sec. 3.4 indeed we have in no case found points correspond¬ 
ing to a mass ordering that would be forbidden if the sum rule held exactly. More generally 
we have seen that the running has in many cases a visible but not a dramatic effect. The 
simple and intuitive reason for this is that the parameters in the neutrino sector are generally 
known to run relatively weakly (although exceptions do exist, see Ref. |46 for an example). 
Thus, even though the sum rules are in reality not anymore valid at the low scale, the running 
effects are sufficiently weak that the sum rules are nevertheless approximately fulfilled for all 
the points displayed, and thus their predictions are not spoiled. The small differences seen are 
negligible compared to the uncertainties coming from nuclear physics, which are however still 


not big enough as to destroy the testability of many groups of neutrino flavour models 115 17 


Hence, we have shown that the RGE effects do not change the qualitative predictions of 
the sum rules, but it should nevertheless not be neglected because they can even have an 
impact for a small mass scale. Especially in the regime with a large mass scale, we have 
shown that the running effects do broaden the allowed region whereas the absence of a visible 
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broadening should be regarded as an exception where “accidental” cancellations take place. 

Thus, while there may be further model-dependent corrections present in case a neutrino 
flavour model yields a sum rule, at least RGE corrections do not change the qualitative 
predictions of the sum rules. In most case, even the quantitative predictions are hardly 
changed, in particular when taking into account that the nuclear physics uncertainties will 
always dominate in a realistic measurement. In turn, the predictions by a certain sum rule are 
safe up to possible model-dependent effects, whose size can however typically be estimated 
or even computed exactly for a given realistic flavour model. 

5 Summary and conclusions 

We have presented the first explicit and systematic study of the effect that radiative correc¬ 
tions have on the validity of neutrino mass sum rules. Since sum rules are able to yield very 
concrete predictions that are realistically testable with near-future experiments, it is impor¬ 
tant to take into account possible modifications if we are to truly put the models developed 
over more than a decade to the test. We have started this endeavour by numerically com¬ 
puting how the regions in the parameter space allowed by certain sum rules are affected if 
renormalisation group running is taken into account. 

After briefly reviewing the most general form of a neutrino mass sum rule and a discus¬ 
sion of the general effect of renormalisation group running, we have explicitly computed the 
resulting allowed regions for all neutrino mass sum rules known if we assume the rules to 
hold exactly only at the seesaw scale, while correction terms appear when going to lower 
energies. The concrete settings we have used were a Standard Model-like scenario (where 
running effects are expected to be very small) and two scenarios corresponding to the min¬ 
imal supersymmetric Standard Model (with tan /3 = 30 and 50, where we expect running 
effects to become stronger with larger tan/3). While we have explicitly verified these general 
tendencies, our results nevertheless show that the predictions derived from neutrino mass sum 
rules, although visibly changed by the corrections, are nevertheless quite stable due to the 
smallness of the effect (this holds unless the running was unusually strong). Three sum rules 
do not run because of cancellations in the /3 functions, or at most by such a small amount that 
the resulting changes in the prediction regions in the m—\m ee \ plane are basically invisible in 
the plots. In fact, not only experiments looking for neutrinoless double beta decay have an 
impact on neutrino mass sum rules. If accelerator experiments determine the neutrino mass 
ordering some cases are directly excluded. Furthermore the LHC can shed light on the ques¬ 
tion whether one has to apply the MSSM or the SM /3 functions for the neutrino parameters. 
This is quite a difference because not only the size of the corrections differs but also the sign 
of the corrections changes. 

Our findings considerably strengthen the position of neutrino flavour models featuring 
mass sum rules, since the predictions derived prove to be relatively insensitive to radiative 
corrections. This leads to a big advantage of such models compared to others not predicting 
any correlation between observables. The only caveat, apart from having a setting where the 
running is very strong, is that some concrete models may induce other big corrections that 
are completely unrelated to the running effects discussed here. While such effects may still 
be able to change the regions predicted by that specific sum rule (or maybe to even entirely 
destroy their validity) in that particular setting, typically both their origin and size would be 
clear in a concrete model, to the point that their strength may even be computed or estimated 
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at least. 

Thus, our results show that the most generic corrections one could think of are, in fact, 
not a problem for neutrino mass sum rules. These types of correlations hence exhibit a strong 
handle that can be used to realistically probe many neutrino flavour models already with 
upcoming experiments on neutrinoless double beta decay, without the need to wait for the 
precision era in neutrino flavour physics. 
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A Taking the square root of a complex number 

There are some subtleties in treating sum rules which include the square root of the masses. 

For a positive real number x one has to take both possibilities of the sign of the square 
root into account, i.e. yfx = ±| y/x\. For a complex number z = p e' x , x £ [0, 27r], in turn, one 
encounters further subtleties. For example, one could either define \fz = \yfp\ e 1X//2 , where 
X £ [0,27r], or one could alternatively define yfz = \yfp\ e 1 */ 2 , where x £ [—'7r, 7r]. In the first 
case the result lies within the upper half of the complex plane, where lm(y/z) > 0 whereas 
in the second example the result is in the right half of the complex plane, with Re(-^A) > 0. 
Depending on the definition, the final results will differ. And special care has to be taken that 
in a numerical setup the two definitions are not messed up. For instance, REAP suggests the 
convention Im.(yfz) > 0 while Mathematica uses Re(-v/F) > 0. 

To cover the whole complex plane one should furthermore consider solutions with a neg¬ 
ative sign of the square root. In our first example, this means that we also have to consider 
y/z = — \y/~p\ e 1 */ 2 . In the following we will employ the definition of the square root according 
to our first example: yfz = Tly 7 /?) e lx / 2 , where x £ [0, 2ir\. Especially in the case of mass 
sum rules which include the square root of the complex neutrino masses a proper definition of 
the square roots is essential since the phases of the neutrino masses have a physical meaning: 
they are the physical Majorana phases. 

An example for a mass sum rule including square roots of the masses is proposed in 
For p = 1 the mass sum rule 10 (see Sec. |4.10[) reads 


12 


y/tin - i \7^3 = 2\/^2, 


(A.l) 


where the masses are all complex. They depend on the complex parameters a, b defined in 
the model 


m\ = (a + b ) 2 , 
m2 = a 2 , 
m 3 = -(a - 6) 2 . 


(A.2) 

(A.3) 

(A.4) 


To get a graphical representation of the sum rule we can, e.g., choose the mass fh 3 to be 
real and positive, fh 3 = m 3 , since we can absorb one phase as a global phase factor. The 
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Figure 3: Graphical representation of the sum rule yjfh\ — i \fm 3 = The small red 

dashed circles represent the left hand side of this equation, the big blue circle represents the 
right hand side. 
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phases of rh\ and m 2 are then the physical Majorana phases. The graphical representation 
is given in Fig. [3j The red dashed circles represent the left hand side of Eq. (A.l), while 


the blue circle represents its right hand side. We have taken into account both possible signs 
for y/rh 3 , which correspond to the centres of the small red circles with radius \\/mi\. The 
big blue circle with radius \ 2 y/ifi 2 \ is centred around the origin. The sum rule is fulfilled if 
and only if the circles have an intersection. If we consider only the positive solution of y/rhz, 
the intersections of the circles in the half-plane where Ke(y/m^) < 0 are absent. Since the 
angles in the triangles formed by the intersections of the circles are related to the Majorana 
phases whose interval is the whole complex plane, we would miss two physical solutions. As 
the circles have four intersection points, we therefore conclude that there are four solutions 
for the Majorana phases which are in accordance with the sum rule, from which only two are 
physical (the other two solutions give the same results). 

From this construction we can as well convince ourselves that the three values of i] from 
sum rule 10 all give the same result. First of all, that the two cases i] = 0 and r] = 2 are 
equivalent is obvious since by construction we have chosen as the center of the red circles 
±\/m 3 anyway. The third case with 77 = 1 can be rewritten to 


\/-rh\ + v^3 = 2\/— 


m2 


(A.5) 


which just mirrors the blue and red circles along the horizontal axis (it adds it to the Majorana 
phases). So this simply interchanges the two physical and the two unphysical (redundant) 
solutions with each other. 

Similar considerations can be done for the other sum rules involving square roots, such 
that there could be equivalent sum rules with additional signs and factors of i. But in this 
study we have quoted only the sum rules we have found in the literature and the underlying 
model fixes the concrete form of the sum rule and hence we do not claim that we list all 
possible mass sum rules. 

Regarding the sum rules which do not include square roots of masses we only obtain two 
solutions for Majorana phases, since we chose m 3 = m 3 to be positive via the redefintion of 
the phases. Hence we only find one circle around m 3 in the right complex half-plane with 
Re(m3) > 0. 

In principle more sum rules can arise when taking different signs of the square roots of the 
masses into account. For all the reviewed sum rules which include square roots of the masses 
we have checked that there is only the quoted combination of signs of the square roots which 
leads to a valid sum rule. This means that there is only one possiblity to form a triangle 
when we interpret the sum rule geometrically. 
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